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Amplitude spectra are deduced for trains of leader pulses and of K pulses from assumed 
models of the field or current variations in the lightning discharge. A comparison is made 
with experimental observations. A size distribution of K pulses relative to return stroke 
pulses is established. Finally a comparison is made between the spectra of atmospherics 
from leader, K and return-stroke pulses. These peak respectively at 20, 8, and 5 kHz, 
with corresponding relative peak amplitudes of 1:2:10. 



1. Introduction 

It is well known that if the electrostatic field 
change due to a close lightning flash is recorded, 
the variation in field with time is comparatively 
smooth [Pierce, 1955a] except for the abrupt change 
associated with the return stroke of the discharge 
to earth. Since the electrostatic field varies with 
M, the electric moment, 1 it might be anticipated 
that radiated electromagnetic fields (proportional to 
d 2 M/dt 2 ) would be of significant amplitude only 
during the return stroke. However, when close 
records are made with equipment having a high 
time-resolving capability or under frequency-response 
conditions that accentuate sudden variations [Ivita- 
gawa and Brook, 1960], it is found that rapid 
changes in M occur during phases of the discharge 
other than the return stroke. The total variation 
in M during these rapid changes is small, but is 
sufficiently fast for d 2 M/dt 2 and, in consequence, 
the radiated fields to be appreciable. It is the pur- 
pose of this paper to consider the VLF pulses gener- 
ated during some of these fast changes. Figure 1 
illustrates their relationship to the overall electro- 
static variation. The upper curves (a) show the 
change of field due to a discharge at 100 km as 
measured with equipment having a frequency 
response of say 1 Hz — 1 kHz; the lower curves (b) 
are for a frequency response of perhaps 1—50 kHz. 
Variations (a) are mostly electrostatic; variations 
(b) dominantly radiation. 

For a cloud discharge, the superimposed pulse 
structure most important at VLF is associated with 
so-called K changes; i.e., small rapid field changes. 
These are generated as a result of the inhomo- 
geneous distribution of charge in the cloud. The K 
pulses are most pronounced in the later stages of a 



1 The electric moment, M, of a .thundercloud is denned as 2 Xqh where q and h 
are the charge and heighi respectively of an element of charge, and the summation 
extends over all the charges associated with the cloud. The quantity M is in 
common usage among physicists and meteorologists. Radio engineers often 
prefer to employ "currenl moment" that is twice (making allowances for the 
image) the product of the currenl flowing and the length it energizes. If the 
tieighl variation is small the current moment is equal to dM/dt. 



cloud discharge; there is evidence that the pulse 
structure immediately following the initiation of a 
cloud flash is somewhat different [Kitagawa and 
Brook, 1960;Ishikawa, 1961]. 

Well-marked K pulses also occur during the inter- 
vals between the return strokes of a flash to earth, 
and following the final return stroke. For the dis- 
charge to ground, the stage preceding the first 
return stroke is usually known as the leader stage. 
It also contains superimposed pulses. Those most 
familiar, which may be termed leader or /. pulses, 
correspond to the luminous steps produced during 
the progress of an initial leader from the base of the 
cloud to the ground [Schonland, 1956]. Pulse 
astructure is also associated with leader development 
from its initiation within the cloud until its emergence 
from the base. There are suggestions that this pulse 
structure does not correspond to true stepping 
[Clarence and Malan, 1957] and that the latter is 
primarily a phenomenon of leader progression 
through free air. 

This paper is concerned with the VLF fields 
associated with the radiation of K and L pulse 
structure. The following approach is adopted. 
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Figure 1. Recordings of a typical cloud flash and multiple 
flash to earth. 
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If the total field variation as a function of time t is 
defined as E(t), then the amplitude spectrum as a 
function of frequency is given by the modulus of the 

J -J- 00 
E(t) exp 

( — iost) -dt and co is the angular frequency. Then 
E(t) may be deduced empirically from field observa- 
tions or derived from estimates of current variations; 
amplitude spectra can then be obtained. These are 
compared with observed amplitude spectra, allow- 
ances being made, where necessary, for any modifi- 
cations due to propagation. Since atmospherics 
are usually recorded at some distance from the flash, 
this implies an extrapolation back to source 
conditions. 

2. Leader Pulses 

A lightning flash is initiated by a leader streamer, 
usually carrying toward the earth negative charge 
from that accumulated near the base of the thunder- 
cloud. It is believed that, having left the cloud, the 
leader maintains an almost continuous advance to- 
ward the ground. However, as regards VLF atmos- 
pherics, the superimposed high-current surges that 
occur, associated with stepping, are more important. 
Observations of these surges, made photographically 
by Schonland [1956], suggest that the charge travels 
about 50 m in less than a microsecond. The bright- 
ness of the channel lasts only a few microseconds, 
but there is an appreciable pause between surges, 
lasting between 40 and 100 jusec. Schonland's 
observations show a mean pause time (r) between 
the luminous steps of about 50 iisec with a standard 
deviation c^23 Msec. A train of leader pulses ob- 
served electromagnetically lasts about 3 msec on the 
average, and contains some 60 pulses. Observa- 
tions made by Kitagawa [1957] and Kitagawa and 
Brook [1960] give 7=70 jusec with o-=50 /xsec. The 
intervals between pulses are observed to decrease as 
the leader approaches the ground [Ishikawa et al., 
1958]. 

If one assumes that the pulses are identical in 
form and that the pause times are normally distrib- 
uted about some mean r with standard deviation c, 
it is possible to derive an expression for the mean 
amplitude frequency spectrum, which we shall de- 
note by |L(w)|. The derivation is as follows: 

If E(t) is the form of a single pulse and F(w) its 
Fourier transform, then the spectrum of a series 
of N pulses, occurring at times T x T 2 . . . T n is 
given by 

S(«) = S exp {-io>T n )F{u) = F(o>) £ exp (-ia>T n ). 

71 = 1 71=1 

The mean power spectrum for such a series is given by 



N 



S(u)S*(a)=F(<*)F*(<*)^ exp (-icT n )J2 exp (icoT n ). 



Expanding the summations under the bar yields 
2V+S exp [iw(T n+1 -T n )]+exp [-ia>(T n+1 -T n )] 



+ 2 exp [iu(T n+2 -T n )]+exp [-iu(T H+i -T % )] 

71=1 

+ further terms 
and, writing T n+ i— T n = T+e n , this reduces to 



N+(N— l)[exp (ioir) exp (icoe w ) 



+ exp (— iayr) exp (— iwe n )] 



+ (A r — 2)[exp (2icx)r) exp io)(e n + e n+ i) 



+ exp (— 2iorr) exp — ia>(e n + e n+1 )] 

+ further terms. 

Assuming the e n all obey the same distribution 
law and are statistically independent, Large [1957] 
has shown that 



exp io)(e 1 + e 2 + . . . e m ) 

° p(e) exp (ioe)del =[P(«)f 

where p(e) is the probability of lying in the range 
e to e+rfe. If AT>>1 and the distribution of e n is 
symmetrical about zero, it follows that 

In actuality, the distribution of e n is skewed with an 
excess of negative values over that expected from 
a standard distribution of the t's, but other approxi- 
mations made in this treatment (e.g., simply the 
fact that we are deriving an average behavior for 
all lightning flashes) render the error in this assump- 
tion negligible. For a normal distribution of the 

r~ ( T 2 o? r \ 
r n , P(o)) = exp — • This gives the mean ampli- 
tude spectrum 

\L(u)\={S(co)S*(co)}i 



Lcosh (|a- 2 co 2 )— COS CO 
= |F(«)|tf*0(«). 



f sinh (jo-V) 7 



The form of 6?(co) is shown in figure 2 for three sets 

of values of 7 and o. It can be seen that when a is 
of about the same order of magnitude as 7 (as is 
experimentally observed), G(u>) does not differ 
greatly from unity, and hence the spectrum |i(co) | of 
the series of pulses is almost proportional to that 
\F(u>) | for the single pulse. However, if a is small, 
and the stepping has therefore a well-defined repeti- 
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Figure 2. The junction (J(w). 



lion frequency, 6?(co) has pronounced maxima at the 
repetition frequency and its harmonics. 

As has already been mentioned, the form of 
F(oS) depends on the form assumed for E(f). Figure 
3 shows two curves derived by the above method 
and compares them with the form (3) given by 
Watt and Maxwell [1957] and with an experimental 
curve derived from the data given in the appendices 
to Steptoe's thesis [1958]. Curve 1, labeled Pierce, 
was derived using Schonland's values for t and a, 
and N i= 8 as suggested by his data. The expression 



n 



E{t)=K[h\t exp {-hlt 2 )+hlt exp (-h 2 2 t 2 

as deduced empirically by Pierce [1955a], was used, 
with h 1 = 3.5X!0 5 sec" 1 and h 2 =l.2X10 5 sec" 1 and, 
since radiation fields are being considered, K is 
inversely proportional to distance for short distances 
from the source: at 100 km K= 2.7X10"' 



v m 



sec-. 



The second curve is derived from the current form 
for a leader pulse suggested by Steptoe, that is, 



J(0 = 2X10 3 [exp (- 



-at)— exp ( 

00^=4X10- 
dt 

exp (— fit)] at rf=100 km. It is assumed, following 

Steptoe, that a step of length 100 m is energized by 



E(t)=-±-~ r 
4x€ c l d 



-j3t)] A, which gives 
[ — a exp (—cd)-\-f} 



the current I (J). Steptoe suggests e*=2X10 4 and 
0=2 X10 5 for the parameters. In this model, the 
current is still appreciably different from zero 50 
jusec after the pulse begins, which seems unreason- 
able in the light of Schonland's data. 
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Figure 3. Spectra at 100 km of a series of leader pulses. 
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Figure 4. Size distribution of K pulses relative to return-stroke pulse. 



The discrepancy of about an order of magnitude 
between the theoretical and the observed results is 
not entirely unexpected. Steptoe's results were 
obtained in conjunction with observations by the 
British CRDF atmospherics locating network [Ocken- 
den, 1947]. With this system, only atmospherics 
larger than a triggering threshold are fixed. In 
addition, Steptoe's equipment responded prefer- 
entially to atmospherics of large amplitude. It is 
believed that these experimental conditions imply a 
bias toward stronger signals. Thus curves 1 and 2 
represent average behavior while curve 4 is for the 
strongest leader disturbances only. 

3. K Pulses 

A similar analysis can be applied to the K pulses. 
The size distribution of these pulses has been 
experimentally studied by Kitagawa [1962]. His 
results can be expressed in terms of the ratio r 
between amplitudes of the K pulses and that of a 
typical return stroke in a flash to earth. The 
following empirical relation applies: 



dP(r) 



1 



dr 



-5/2 



8V2tt 



exp 



< 



\ 8r 






-F(r), 



where dP(r) is the probability of values between r 
and r+dr. This distribution is plotted on figure 4; 
it gives 0.05 as the most probable value of r and 
0.25 as the mean value. This agrees fairly well with 
the values 0.05 <r< 0.2 given by Kitagawa and 
Brook [1960] and r=0.1 given by Ishikawa [1961]. 

The form of a if pulse is not well known, but 
assuming a double exponential current surge i(t)=i 
[exp (— at) — exp (—fit)] as suggested by Steptoe 
[1958] and making the modification that the velocity 
of advance of the surge along the channel is given 
by v(f)=v exp (—yt) gives 



E{ty- 



__1 1 2^0 

47T€o C 2 d 7 

+ (a+7)exp- 



-a exp (— at)-\-(3 exp (— pt) 
■(a+y)t-(P+y) exp-(fi+y)t]. 



Combining the data of Steptoe and Ishikawa [1961], 
we take i =16,000 A; ^ =2X10 7 m/sec, a=5X10 4 
sec -1 , 0=2X1O 5 sec" 1 , and T =4X10 4 sec" 1 . These 
parameters correspond to a peak current of about 
7500 A reached in a time of 9 /xsec; the channel 
length at the instant of peak current is approximately 
150 m, and the asymptotic limit is 500 in. 

It is found [Kitagawa and Brook, 1960] that 
r=8.5 msec with a=6 msec for K pulses. The 
duration of a series of K pulses depends on the type 
of flash. The junction stage, between return strokes 
of a flash to earth, generally lasts about 60 msec, 
while the final stage of such a discharge will continue 
for about 150 msec. A cloud discharge usually lasts 
about 250 msec. Thus a series of K pulses may 
contain anything from 6 to 30 pulses. Since cloud 
discharges are most common, we shall take N=25, 
giving N* = 5> for this analysis. 

The semitheoretical curve thus obtained is plotted 
on figure 5 and compared with an experimental 
curve deduced from the data in Steptoe's thesis. 
Again, the bias due to triggering makes the experi- 
mentally measured amplitudes considerably larger 
than those given by the theoretical curve. 
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Figure 5. Spectra at 100 km of a series of K pulses. 
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4. Discussion 

Amplitude spectra have been deduced for a series 
of K pulses and for the train of L pulses associated 
with a stepped leader. In order to relate these to 
the size of a return-stroke pulse, figure 6 has been 
prepared. The return-stroke data are obtained from 
a companion paper [Dennis and Pierce, 1964]. 

It can be seen that the spectra for a first return- 
stroke pulse, a series of K pulses, and a series of L 
pulses have maximum amplitudes in the ratio 
10:2:1, occurring at 5, 8, and 20 kHz, respectively. 
Thus, even in the VLF range, the contributions to 
radio noise from processes other than return strokes 
are appreciable. This is particularly so for tropical 
areas, where cloud discharges are perhaps 10 times 
as frequent as flashes to earth. At high frequency, 
of course, it is well established [Clarke, 1960] thai 
the return stroke is of little significance in the 
production of radio noise. Previous work [Pierce, 
1955b, 1957; Watt and Maxwell, 1957] has pointed 



out that the leader can produce appreciable noise 
at VLF and LF. The present paper shows that the 
K pulses are considerably more important that the 
leader disturbance. Within the VLF band, figure 6 
indicates that upon the average the signals go in 
ascending order of magnitude as leader, K change, 
return stroke. However, the dynamic range at the 
source is large for all three types of disturbance so 
that many K changes and even a few leaders are 
greater than some return strokes. Since K changes 
occur much more frequently than return strokes their 
importance as regards VLF noise is evident. 



We are indebted to Dr. B. J. Step toe for being 
allowed to use the data contained in his thesis. 

This paper summarizes work described in greater 
detail in a report by Pierce, Arnold, and Dennis 
[1962] and supported under U.S. Air Force Contract 
AF 33(657)-7009. We are grateful to the Air Force 
authorities for permission to publish this paper.J , &J 
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Figure 6. Comparison of spectra at 100 km due to various processes in a lightning flash. 
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